Precision Measurement of the ^He Half-Life and the Weak Axial Current in Nuclei 
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Studies of ®He beta decay along with tritium can play an important role in testing ab-imtio 
nuclear wave-function calculations and may allow for fixing low-energy constants in effective field 
theories. Here, we present an improved determination of the ®He half-life to a relative precision 
of 3 X 10"*. Our value of 806.89 ± diktat 

-0.19 Byat ms resolves a major discrepancy between 
previous measurements. Calculating the statistical rate function we determined the /t-value to 
be 803.04tg;2^ s. The extracted Gamow- Teller matrix element agrees within a few percent with 
ab-initio calculations. 
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Precision measurements of electroweak processes in 
liglit nuclei can provide important tests of our un- 
derstanding of electroweaic interactions in the nuclear 
medium. Many interesting problems - ranging from solar 
fusion to neutrino interactions and muon and pion cap- 
ture processes - depend on their correct modeling and 
calculation Recent progress in numerical techniques 
enable precise, ab-initio calculations of wave functions for 
light nuclei to be performed with the nucleon-nucleon in- 
teraction and without assuming a frozen core of inactive 
particles 

The allowed weak nuclear decays driven by the ax- 
ial current - called Gamow- Teller decays - have histor- 
ically played an important role in testing wave func- 
tions because the main operator has a simple spin and 
isospin structure and does not possess any radial com- 
ponent. Systematic comparisons with shell-model wave 
functions showed that in order to reproduce observations 
tlic value for the weak axial coupling constant, gA, had 
to be 'quenched'. For the sd-shell nuclei this difference 
amounted to about 30% with respect to that measured 
in free neutron decay [5, JiJ. In addition, when charge- 
exchange reactions were used to explore a large fraction 
of the Gamow- Teller strength sum rule, evidence also 
showed up for 'quenching of the Gamow- Teller strength' 
0, S|. However, the origin of the quenching is not com- 
pletely clear. Reference [9] and others Q have shown 
that, as shell model calculations are allowed to intro- 
duce higher and higher excitations, the need to renor- 
malize operators disappears. But it has also been pointed 
out that meson-exchange currents (mediating, for exam- 
ple, nucleon-delta excitations) could be responsible for at 
least some of the apparent quenching of strength . 

The decays of '^H and of ^He are special because these 
systems are light enough that the corresponding ab-initio 
calculations can be performed with precision. In particu- 
lar, Refs. 0, [m and later 12, have shown that, using 
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FIG. 1. (Color online) Compilation of the measured ®He half- 
lives found in literature [l5l- [2lj |. The dashed, blue band shows 
the half-life adopted in Ref^~[2^ from the average of the two 
values found in Ref. [l^ . [20| and used in compilations ever 
since. The inset shows the five values with uncertainties below 
1% [l^ - [20l ] with the red band depicting the value for the ^He 
half-life obtained in this paper. 



the case of '^H to fix nucleon-delta excitations, they can 
reproduce the ft-vahie for ®He to within a few percent. 
These two decays, then, can play an important role in 
testing the accuracy of nuclear wave-functions calcula- 
tions [J, [nl [13, or as suggested in Ref. 13], in fixing 
low-energy constants in effective-field theory calculations 

In this paper we present a high-precision experimental 
determination of the half-life and ft- value for ^He. Ex- 
cept for a small branch of ~10~^ [1J| the beta decay of 
^He proceeds exclusively to the ground state of ^Li with 
an endpoint of 3.5 MeV. Its half-life has been previously 
determined by several works compiled in Fig. [TJ As can 
be seen, the values spread over a range much wider than 
expected from the claimed uncertainties which makes the 
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currently reported average and precision of 806. 7 ±1.5 ms 
[2^ unreliable. Averaging the five values shown in the in- 
set in Fig. [T] with uncertainties below 1% and scaling the 
uncertainty by the square root of the x^/dof - as advised 
by the Particle Data Group in such cases - results 
in 800.6 ±2.0 ms. As emphasized by the author of the 
last precision measurement in 1982 '20| this needs to be 
resolved by improved and higher precision experiments. 

We produce ^He using the tandem Van de Graaff accel- 
erator available at the Center for Experimental Nuclear 
Physics and Astrophysics of the University of Washing- 
ton via the reaction ^Li(^H,^He)^He. A detailed descrip- 
tion of the ^He source, which can deliver ~10^ atoms/s 



to experiments, can be found in Ref. [2^. The °He 
atoms are transferred from the source through a turbo- 
molecular pump to a low-background experimental area. 
During a period of 8 s we direct the deuteron beam onto 
the lithium target and the outlet of the turbo-molecular 
pump is connected to our measuring volume. For the fol- 
lowing 16 s, thereby completing one cycle, we deflect the 
deuteron beam at the low energy end of the accelerator, 
close a spring loaded valve in front of the measuring vol- 
ume and measure the ^He half-life directly by observing 
the decay curve. 

The measuring volume consists of a 35 mm diameter, 
381 mm long tube made of stainless steel sealed on one 
side by a 254 iim thin copper foil and on the other by 
a spring loaded, viton 0-ring sealed valve. Randomly 
distributed over our measurement period of five days and 
taking up about half of the available data, we inserted a 
19 mm diameter, 283 mm long stainless steel cylinder 
suspended in the center of the tube by four bolts. With 
the stainless steel insert we increased the wall collision 
frequency of ^He atoms by about 80% and used that data 
to check for possible diffusion of the ^He atoms into the 
stainless steel surfaces. 

Directly in front of the copper foil we placed two iden- 
tical, 2.5 mm thick plastic scintillators registering the 
betas from the decay of ^He. The plastic scintillators 
were coupled via lightguides to photomultiplier tubes. 
Their output signals passed through timing filter am- 
plifiers (resulting in ~100 ns long pulses) and discrim- 
inators set to the tail of the electronic noise (cutting 
away ~1% of the electron spectrum). We formed the 
coincidence of those two signals resulting in a single 
25 ns long logic pulse and a total intrinsic deadtime of 
~130 ns, as measured on an oscilloscope. We passed this 
pulse through four gate generators providing signals of 
fixed, non-extendable deadtimes of lengths 1.9819(81) /iS, 
3.9990(81) us, 6.0026(83) ^s, and 7.9758(83) fis. The 
deadtimes were determined following the ^He run using 
the source and pulser method with the pulser sending in 
logic pulses at the stage of the long gate generators (2^ . 

The four gate signals were fed into a CAMAC based 
scaler together with the original coincidence signal and 
the signals from a 1 and 100 kHz clock. The signal from 
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FIG. 2. ®He decay curve for the data with initial rates 
<40 kHz and the stainless steel insert out with the corre- 
sponding residuals. The x^/dof of the fit amounts to a typi- 
cally observed 1578.2/1562. 



the 1 kHz clock also triggered the read-out of the scaler 
module via the software package JAM [2^ thereby pro- 
viding 1 ms time stamps to our data stream. 

We took data at various deuteron beam intensities 
to study decay-rate dependent effects. The data were 
grouped into five different initial rate classes: < 40 kHz, 
40 - 50 kHz, 50 - 60 kHz, 60 - 70 kHz, 70 - 80 kHz. 
We corrected the deadtime losses after rebinning on a 
cycle-by-cycle basis by calculating the true rate Rq in 
each time bin from the measured rate R using the mea- 
sured deadtimes Td- Ro = R/ {^-Rrd) [27|. Figured] 
shows the decay curve for the data with initial rates 
<40 kHz and the stainless steel insert out. Also shown 
is a fit and the corresponding residuals with the function 
R{t) = N{exp{-{t-to)/T) + b), where the lifetime t 
and the background b are free parameters while N is set 
by the total number of counts. The fit was performed 
using the modified method outlined in Ref. |28!|. We 
cross-checked our fit results performing both regular 
and maximum likelihood fits. The effects of data treat- 
ment and fitting method were studied using simulated 
data without noticing a particular bias. 

As the systematic uncertainty due to the deadtime cor- 
rection grows rapidly with rate we delayed the starting 
point of our fit in each of the rate groups such that for the 
highest rate in that group the initial rate lies at 32 kHz. 
Table U gives the final results from the averages of the 
different rate groups. The individual deadtime channels 
were combined to yield the average values for the two 
cases of the stainless steel insert in and out of the mea- 
suring volume. These two values are used below for an 
estimation of the potential diffusion of ®He atoms into 
the stainless steel surfaces. 
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TABLE I. List of the difFerent half-lives obtained with various 
deadtimes and for the cases of the stainless steel insert in and 
out. The systematic shifts and uncertainties from Table [TI] 
are not included. 



Insert 




Results [ms] 


Average [ms] 






807.01 ±0.11 




out 


-^4 /is 


807.03 ±0.11 


807.03 ±0.11 




-^6 /is 


807.08 ±0.11 






-^8 /is 


807.02 ±0.11 






-^2 fis 


807.20 ±0.12 




in 


--4 /is 


807.20 ± 0.12 


807.21 ±0.12 




-^6 /is 


807.21 ±0.12 






~8 /is 


807.23 ±0.12 





TABLE II. List of systematic shifts and uncertainties. We 
added the errors in quadrature to obtain the total error. 
Where a second value is given it corresponds to the mea- 
surements with the stainless steel insert in. 



Source 



Shift [ms] Uncertainty [ms] 



Deadtime correction 




0.037 


''He Diffusion 





+0.12/+0.22 
-0 


Gain shift 


-0.19 


0.19 


®Li contamination 





+0 

-0.007 


Background 


0.046 


0.004 


Data correction 





0.01 


Deadtime drift 





0.009 


Afterpulsing 





0.003 


Pile-up 





< 0.005 


Clock accuracy 


0.006 


0.011 


Total 


-0.14 


+0.23 / +0.29 
-0.19 / -0.19 



In Table |TT] we give our estimates of identified system- 
atic shifts and uncertainties. For brevity we will only 
describe the most important ones in the following para- 
graphs. A detailed description will appear elsewhere [29| . 

The effect of the deadtime correction is large. The 
correction shifts the half-life values by —9 ms, —18 ms, 
—27 ms, and —35 ms for the four different deadtime chan- 
nels, respectively. However, the agreement between the 
four values after the deadtime correction lends confidence 
to its validity. The uncertainty on the measured dead- 
times of ~8 ns translates directly into a systematic un- 
certainty on the half-life with its dependence extracted 
from the data. 

Using a helium leak detector we studied the amount 
and time evolution of the diffusion of helium through the 
walls of the measuring volume and the viton 0-ring of 
the valve at its end. The observed amount of diffusion 
through the 0-ring leads to a negligible shift at our level 
of precision. However, there most probably occurs some 
amount of diffusion into the walls. Any time constant 



I/tzoss associated with this loss channel scales linearly 
with the wall collision frequency. From the absence of 
any significant difference between the two results listed in 
Table U] we conclude that tioss 3> Tepjo and the difference 
between the two results amounts to A(1/t) — l/rin — 
1/Tout (-2.8 ± 2.5) X 10-^ ms-i = O.S/noss- We set 
the corresponding Gaussian probab ility density function 
to zero in the non-physical region [23| and calculate an 
upper limit on 1/tioss at 68% C.L. of 2x10"'' ms"^. This 
translates into systematic uncertainties for the insert in 



and out data of to'^^ 



and 



+0.12 
-0 



ms, respectively. 



Examining our highest rate data we have identified 
traces of a small, rate-dependent shift, which is not fully 
accounted for by our deadtime correction. The half-lives 
are shifted towards higher values and we attribute this 
to a negative gain shift with increased rate in the pho- 
tomultiplier tubes. In order to investigate its infiuence 
on our data constrained to 32 kHz, we added a parame- 
ter k to our fits in order to model the effect of a linear, 
rate-dependent shift by substituting in our fitting func- 
tion R(t) -)> R(t) (1 - kR{t)). The resulting shift in the 
half-life due to including the parameter k amounts to 
—0.19 ±0.19 ms. The fact that this parameter is com- 
patible with no remaining rate-dependent effects was also 
observed when plotting the results of our half-life fits as 
a function of starting point of the fit window. The varia- 
tions of the results were in agreement with those expected 
from the loss of statistics. 

A potential contaminant in our system - apart from 
tritium, which will not influence our measurement - is 
the beta emitter ^Li (T1/2 = 838.40(36) ms [s^l with end- 
point ~16 MeV) produced by the reaction ^Li(^H,^H)^Li. 
The lithium atoms are not expected to reach our count- 
ing station but rather get trapped in the lithium target 
or in the walls during the many collisions (^10^) that 
occur before atoms can reach the detection area. Nev- 
ertheless, in separate measurements using both our two 
thin detectors used for the lifetime measurements but 
also a thick scintillator to measure the full energy of the 
betas we scanned the deuteron beam energy below the re- 
action threshold for ^He production via ^Li(^H,'^He)^He 
of 5.8 MeV but above the one for ^Li of 0.25 MeV 
While we still observed production of a beta emitter both 
the energy spectrum with an endpoint of 3.5 MeV and 
the extracted half-life of 810 ±5 ms clearly identifies it as 
^He. The observed production agrees with a rough esti- 
mate of it being produced by the reaction ^Li(n,^H)^He 
where the neutrons stem from ^Li(^H,n). Integrating 
the background-subtracted beta energy spectrum above 
the ^He endpoint and operating at our nominal deuteron 
beam energy we set a limit of 2 x lO"** at 68% C.L. on 
the fraction of possible ^Li contamination. 

Throughout the data taking we performed several 
background runs in which we kept the valve in front of 
the measuring volume closed but otherwise operated the 
experiment just like for the other runs. While we ini- 
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tially saw a significant contribution stemming from ^He 
betas penetrating through the thin-walled stainless steel 
bellows of our roughing pump, we were able to greatly re- 
duce that background by shielding the bellows with lead. 
Combining all the background run data, we still observe a 
small decay structure with a half-life of 507(27) ms and 
an amplitude of 6.3(3) times the value of the constant 
background of 0.8 ± 0.1 Hz. While this is most proba- 
bly still coming from ^He that is being pumped away it 
could also be the result of some beam related activation. 
Regardless of its origin, this time-dependent background 
results in a systematic shift of 0.046 ms with an uncer- 
tainty of 0.004 ms. 

As the result of our measurements with the stainless 
steel insert is dominated by the systematic uncertainty 
due to a potential diffusion of the ^He atoms into the 
surfaces we do not average the two values given in Ta- 
ble m We report the data from our measurements with- 
out the insert as our final result yielding a ^He half- 
life of 806.89 ± O.llstat to^igsyst ms. From this, we pro- 
ceed to determine the ft-valne for the beta decay of ^He 
and extract the corresponding Gamow- Teller matrix el- 
ement. We calculated the Q- value of the decay to be 
3.505208(53) MeV/c^ using the recent '^He mass deter- 
mination obtained in a Penning trap [si'l and the value 
for ^Li [sl]. This corresponds to a 4ct shift compared to 
previously reported values [3J] . The relation between the 
/i- value and the Gamow- Teller matrix element Mqt is 
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(1) 



following the definitions and notation of Ref. [35| . We set 
the parameters ^ns and 5c to zero - or equivalently ab- 
sorb them into the definition of Mqt - and calculate the 
radiative correction to be 1.0365(13)%. We adopted 
the value for the parameters iiT/ (Gy(l -I- A)^)) = 
6143.62 ± 1.66 s from the world av erag e of superal- 
lowed O"*" 0+ nuclear beta decays [35|. The statis- 
tical rate function is given by /* = / F{Z, E)pE{E — 
Eoffi{E)dE = /(I + 5s) where / is the value of the 
integral in the absence of the shape-correction function 
fi{E) and Sg the correction to it when including fi{E). 
Here F{Z, E) is the Fermi function, p and E the elec- 
tron momentum and energy, and Eq the end-point en- 
ergy. We obtain / = 995.224(68) yielding an /i-value of 
803.04to'23 s, where we added the statistical and system- 
atic errors in quadrature. In order to take into account 
the shape correction we performed shell model calcula- 
tions using the Cohen-Kurath interaction [36] and with 
the PWBT interactions [37]] adjusted to either reproduce 
the experimental Gamow- Teller matrix element or the 
weak magnetism term, which in Holstein's notation 38 1 
is 6 = 68.4(7), determined from the width of the 0+ — 1+ 
transition in ^Li [s^. Both adjustments result in almost 
identical terms for the statistical rate function and we 
obtain /* = 997.12(58). From this we calculate the ex- 



perimental value for the Gamow- Teller matrix element 
in ^He beta decay as |Mgt|= 2.7491(10)/|gA|- Using 
gA — —1.2701(25) [3] determined from the decay of the 
free neutron, we get |Mgt|= 2.1645(43). 

Given the precision obtained in the measured half-life 
it is worthwhile to carefully include small contributions 
affecting the calculated /t- value. The comparisons so far 
[1, [TTI - [ia | used an /t- value which dates back to Ref. [i^l 
and a statistical rate function obtained from tabulated 



values in Ref. [4l|. By happenstance our half-life mea- 
surement agrees with the value adopted in [2^ . The cal- 
culation of the statistical rate function and radiative cor- 
rection together with the large shift of the Q-value due 
to the recent ^He mass measurements results in an ex- 
perimental Gamow- Teller matrix element also in agree- 
ment with the one used in comparisons so far. How- 
ever, the extraction of the experimental Gamow- Teller 
matrix element now stands on much more solid footing. 
Neglecting the error due to qa the difference between 
the experimental and most recent theoretical matrix el- 
ement [l^l amounts to 0.034(5) with theory overpredict- 
ing the experimental value by 1.5%. The other calcu- 
lations are within -5% [llil|, El- Only [3| gives the 
calculated width of the analogous Ml transition in excel- 
lent agreement to within the experimental uncertainty of 
2%. Interesting, accompanying measurements to our re- 
sult would be improved determinations of the Ml width 
and of the muon capture rate on ^Li [4^ to a preci- 
sion of <1%. The latter would directly test a possible 
momentum-transfer dependence of axial currents in nu- 
clei 0. 

In summary, we have performed the most pre- 
cise measurement of the ^He half-life of 806.89 ± 
O.llstat ^019 syst ms thereby iniproving the precision over 



the currently reported value [22] by a factor of 6. Our 
result is in good agreement with two of the previous five 
values [ll, 113 with precisions of less than 1% but devi- 
ates from the three others by up to 8.6cr [3, 17, \\^ . 
Because the possibility of diffusion out of the target 
was not directly addressed we speculate that this may 
be the cause of the discrepancy. Calculating the sta- 
tistical rate function we determined the /t-value to be 
803.04to:23 s. The extracted Gamow-Teller matrix ele- 
ment of ]Mgt1= 2.1645(43) agrees within a few percent 
with ab-initio calculations using the weak axial coupling 
constant gA measured in free neutron decay. Our precise 
determination allows for improved comparisons between 
theory and experiment and may allow using ^He in addi- 
tion to ■^H to fix low-energy constants in the effective-field 
theory description of the electroweak processes. 
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